Droplets in Disordered Metallic Quantum Critical Systems 



O 

o 

(N 
00 



X3 
O 

o 



> 
o 



o 



X 



A. H. Castro Neto'"'* and B. A. Jones 
"'Department of Physics, Boston University, Boston, MA 02215 
IBM Almaden Research Center, San Jose, CA 95120-6099 
(February 1, 2008) 

We present a field theory for a structurally disordered magnetic system coupled to a metallic 
environment near a quantum critical point. We show that close to the magnetic quantum critical 
point droplets are formed due to the disorder and undergo dissipative quantum dynamics. We show 
that the problem has a characteristic energy scale, the droplet Kondo temperature, that determines 
the crossover energy scale from weak to strong coupling. Our results have direct significance for the 
Griffiths-McCoy singularities of itinerant magnets. 
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The behavior of itinerant magnets close to quantum 
critical points (QCP) has been a subject of intense re- 
search in the last few years. It has been found experimen- 
tally that a broad class of systems show anomalous metal- 
lic behavior in the paramagnetic phase Q. Since a large 
number of such systems are structurally disordered due 
to chemical substitution the question arises of the impor- 
tance of disorder for the understanding of the anomalies 
observed in the experiments. We have proposed recently 
that the anomalous behavior observed in some of these 
systems can be understood in terms of Griffiths-McCoy 
singularities close to a QCP These singularities oc- 

cur in the context of percolation theory on a discrete lat- 
tice when clusters of spins tunnel quantum mechanically. 
In the presence of a metallic environment we have shown 
that electrons scatter against clusters leading to a clus- 
ter Kondo temperature associated with their dissipative 
quantum dynamics We found that while dissipation 
freezes the clusters when the system is close enough to 
the QCP, there is still possibility of quantum behavior in 
a large region in the parameter space around the QCP. In 
this case small Kondo temperatures can be obtained even 
for relatively small clusters dropping the requirement of 
cluster "rarity" as a condition for anomalous magnetic 
behavior These results put our approach in proxim- 
ity to the single ion Kondo disorder theories Q. 

In a recent paper Millis, Morr and Schmalian (MMS) 
proposed that a single isolated local perturbation close 
to a QCP of an itinerant Ising magnet produces large 
droplets with internal structure that are blocked from 
tunneling due to dissipative effects. While that theory 
has similar features to the Grifhths-McCoy scenario pro- 
posed by us, their theory is supposed to be valid for clean 
critical systems with a vanishing small amount of local 
defects. Experimentally this scenario can be realized by 
introducing a very small amount of impurities into clean 
stoichiometric compounds like CeRu2Ge2 [D or CePd2Si2 
and driving the system close to a QCP by application 
of hydrostatic pressure. 

In this work we study the problem of a disordered alloy 
such as Ui-x^x^^z [pi or UCu5_a;Pdj, with a finite 



density of defects or impurities and where the distance 
from the QCP is controlled by chemical substitution. We 
show that, for the same field theory studied by MMS, a 
finite density of impurities leads to Griffiths-McCoy sin- 
gularities in a region close to the QCP. We show that in 
the presence of disorder the droplet is structureless and 
its average size is set by the magnetic correlation length, 
^. More importantly, we demonstrate that the droplets 
have a finite droplet Kondo temperature, Tr- , that varies 
continuously with the distance from the QCP and van- 
ishes when the system is sufficiently close to it. These 
results are in agreement with our previous microscopic 
analysis of the disordered Kondo lattice model Q but 
also apply to other models of strongly correlated elec- 
trons such as the Hubbard model [ p^ . 

The starting point of our analysis of droplet formation 
close to a QCP is the Hertz action for a critical itiner- 
ant magnetic system in d spatial dimensions ||l0|] (we use 
units such that h = Ub = 



n,q 



r) |v?(w„,q)|^ 



where aj„ and q are the Matsubara frequency and mo- 
mentum, respectively, 7 is the coupling between the or- 
der parameter (/^(x, r) and the particle-hole continuum, 
Tq gives the momentum dependence of the dissipative 
coupling {Tq — q'~^ with ^ = for antiferromagnets, 
C = 1 for clean ferromagnets and C, = 2 for disordered 
ferromagnets) . The distance from the QCP is controlled 
by r. In the ordered phase r < indicating an instability 
towards long range order and at the QCP we have r = 0. 
In this work we focus entirely in the paramagnetic phase 
where r > so that we can parametrize r = Disor- 
der is introduced into the problem as a random variation 
of r in real space: 







dr (5r(x) 'p^(x, ' 



which is assumed to be Gaussian distributed with width 
M, that is, the probability distribution is given by: 



1 



P{5r) oc exp 



1 

4u 



dx ((5r(x))^ 



(1) 



so that (5r(x) = and 



(5r(x)(5r(y) = u (5'^(x — y) 

where the average is calculated with (0) . The reason for 
the appearence of droplets in the problem is that the 
Gaussian distribution (0) allows for local values of (5r(x) 
(the tails of the distribution) such that r — (5r(x) < 0, 
that is, it allows for local order even in the absence of 
long range order. These locally ordered regions in a sur- 
rounding paramagnetic media are called droplets. 

Since we are not interested in one particular realiza- 
tion of the disorder we study the average free energy us- 
ing replicas We introduce n replicas of the order 
parameter Lpa with a — 1, ...,n and calculate the average 
free energy, F = (Z" — taking the limit of n — > 
at the end of the calculation. It can be easily shown that 
a new term is generated in the quantum action: 
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dx 
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Notice that the disorder not only generates interactions 
between fields in different replicas but also couples the 
fields in the imaginary time direction. Our calculations 
can also be carried out with a non-linear term of the 
form g(/3'*(x, r) with no fundamental change in the results. 
Thus, in order to keep the discussion simple we have 
dropped this term. 

We first study the problem of droplet formation by 
investigating the static, classical, part of the action. In 
this case ipa (x, r) — ipa (x) where tpa is obtained from the 
variational solution of the static part of the action: 

-VVa(x)+r^a(x)-u/?7Aa(x)^1/^b'W =0. (2) 

b 

This equation is quite revealing. If the problem were 
to be replica symmetric, that is, V'a(x) = V'o(x) for all 
values of a, the last term in (||) would scale like nilj^{'x.) 
and would vanish in the limit of n ^ 0. This would imply 
that the only solution in the paramagnetic phase (r > 0) 
is the trivial solution -00 = 0. Thus, in order for droplets 
to form in the paramagnetic phase one needs the replica 
symmetry to be broken. Here we follow Dotsenko ^ 
who studied the classical problem in detail and assume a 
non-trivial replica solution such that for a = l,...,k we 
have V'a(x) — '0fc(x) while ?Aa(x) = for a = k + 1, ...,n. 
Here A; > 1 is an integer that determines the degree of 
the symmetry breaking process. Using this particular 
solution we see that (0) can be rewritten as: 



This equation is non-linear Schrodinger equation and can 
be thought as the equation for a classical particle moving 
in the potential shown in Fig.^. 
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FIG. 1. Replica potential. 

A naive conclusion from this discussion would be that 
the "energy" would be minimized for ip = ±oo, a solution 
that is clearly unphysical. It should be kept in mind 
that the replicated action is not the actual free energy 
of the problem; the latter is only obtained after we take 
the limit of rt ^ at the end of the calculation (this 
fact is a standard consequence of the application of the 
replica method [0). Dotsenko has shown that the stable 
static solutions of this problem are the maxima of Fig||. 
Mathematically this can be proved by showing that all 
the eigenvalues of the fluctuation matrix (the Hessian) 
are real and positive jlj]. The maxima are associated 
with the two possible configurations of the droplet (spin 
configurations pointing up and down). Equation (^) can 
be rescaled if we define 



uf3k 



(f>{Vrx) 



so that 0(2:) obeys a scale independent equation: 
- Wl(j){z) + (I) ~ <j>^ = 



(4) 



(5) 



where z = y^x. The proper boundary conditions are 
(/)(0) — constant and 4>{z — » ±00) — 0. Equation (||) 
has exponentially decaying solutions for 2; ^ 1/ y/r and 
is smooth for x < l/\/r. This can be contrasted with 
the solution found by MMS for a local defect in which 
the droplet is such that for a; < 1/ y/r the droplet profile 
decays like 1/x. It is the 1/x decay that makes the MMS 
problem special . Here such a behavior does not occur 
in d < 4. Moreover, it is clear that the size of the droplet 
is the magnetic correlation length: 



VVfe(x) + rV'fc(x) - upkiPli^) = . 



(3) 



Thus droplets become arbitrarily large close to the QCP. 

The action for the static droplet can be calculated by 
substitution of (|^) into the original action, and the free 
energy, after summing over the replicas, reads p^ : 
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where — J d'^z(p'^^{z) and V is the volume of the 
system. As pointed out by Dotsenko the non-analytic, 
non-perturbative, dependence of the free energy on the 
disorder strength u for d < 4 shows that the static field 
theory reproduces the classical Griffiths result in a perco- 
lating lattice . Obviously in order to study quantum 
or Griffiths-McCoy singularities we have to allow the 
droplet to tunnel between the two maxima of the poten- 
tial in Fig.0. 

In order to study the tunneling of the droplet we as- 
sume a rigid droplet approximation : 



(^(x,r) = ■(/'fc(x)X(r) 



(6) 



where all the dynamics is encapsulated in X{t). This 
choice assumes that the droplet tunnels as a whole. Di- 
rect substitution of (||) into the action leads to: 
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is the replica free energy for X{t). The coefficients in (Q) 



are: 
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(8) 



that can be associated with the "particle" mass and 
dissipation coefficient, respectively. Here Gq = 
J d'^qq~''\(j)q\'^, 4)q is the Fourier transform of 4){z) and 
periodic boundary conditions are assumed: X{t + j3) = 
X{t). Observe that T describes a Caldeira-Leggett ac- 
tion jist for the motion of a dissipative particle in a 
non-local, non-linear, potential which, for slowly varying 
configurations of X(t) reduces to the potential of Figj^. 
Assuming that most of the time the field is in the equi- 
librium configuration {X{t) = ±1) the problem becomes 
equivalent to the problem of a two-level system coupled 
to a dissipative environment p^ ]. 

The amount of damping in the tunneling of the droplet 
can be estimated by comparing the parameters in (|^). 
We have weak damping when ry ^ loqM while for strong 
damping one has rj ^ ujqM {ujq = 2 is the undamped 



frequency of motion in our dimensionless units). Notice 
that the crossover from weak to strong dissipation occurs 
at r = Tc where (using (ph) 



\2Ei 



2/C 



which indicates that close to the QCP (r < Tc) the 
droplet motion is highly damped. The tunneling split- 
ting between the two configurations A = ±1 is given by 



-SMuial 



where cud is the classical frequency of oscillation. On the 
one hand in the weakly dissipative regime (r > Vc) one 
has uJci ~ Wo and therefore the tunneling splitting is: 



(9) 



where a = I6E1/E2 is a universal constant. On the other 
hand in the strongly dissipative regime (r < Vc) we have 
LOci ~ Muj'^/rj and therefore 



Ar<r,(r) 



4//2 



(10) 



where b — ?>2Ei/ {jGc} is a non-universal constant. For 
antiferromagnets (C — 0) the tunneling splitting is con- 
stant close to the QCP while it vanishes in the case of 
ferromagnets (C ~ 1,2). Another important parameter 
is the Caldeira-Leggett dissipative coupling that is given 
by: 



a{r) = — = (—] 
where we used (^) and defined 

l/(l+C/2) 
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(11) 



(12) 



Observe that the dissipative coupling diverges at the 
QCP. As is well-known this problem has a characteris- 
tic crossover energy scale that can be associated with the 
Kondo temperature, Tr-, of an anisotropic Kondo impu- 
rity model Q . This energy scale separates the region of 
strong and weak coupling and for a < 1 is given by: 



TK{r) « W 



A(r) 
W 



(13) 



where is a cut-off energy scale. For a > 1 we have 
Tk = and the droplet is frozen. Notice that accord- 
ing to ( prj ) the freezing of the droplet occurs for r < rg. 
This indicates that close to the QCP droplets are frozen. 



Observe that vq given in (12) is a non- universal quantity 
which depends on parameters that cannot be obtained in 
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a continuum field theory. Depending on the microscopic 
parameters the freezing of the droplets can occur in the 
region of weak damping if tq > Tc or strong damping 
if ro < Tc, affecting the value of the tunneling splitting 
given in (|^) and (^0|). Thus Txir) is finite and a continu- 
ous function of the distance from the QCP in agreement 
with our previous analysis [|j . This should be contrasted 
with the result obtained by MMS g where Tjf = in 
all the parameter space. In fact for r ^ ro the droplets 
become free to tunnel The phase diagram is shown 
in Fig.| 
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Frozen Droplets Tunneling Droplets 
FIG. 2. Phase diagram as a function of the parameter r: 
the system is ordered for r < 0; droplets are frozen in the 
paramagnetic phase for < r < ro and can quantum tunnel 
for r > ro. 

Direct comparison of ([Tsl ) with our previous results in 
Rcf. m indicates that the two problems map into each 
other if oc where N is the number of spins in the 
cluster. Another interesting consequence of our calcula- 
tion is the sharp contrast between the case of ferromag- 
netic and antiferromagnetic droplets. For antifcrromag- 
netic droples (C — 0) the dissipation coefficient in ( |ll| ) 
scales with while in the case of a clean ferromagnet 
(C = 1) we find a (X indicating stronger damping. For 
a ferromagnetic system with diffusive electrons (C = 2) 
damping becomes even stronger with a oc This indi- 
cates, in agreement with our previous analysis that 
dissipation is more important in the case of ferromagnetic 
droplets than in the case of antiferromagnetic ones. Since 
the great majority of systems studied experimentally are 
of the antiferromagnetic type our results indicate that 
dissipation does not freeze the droplets in a region of the 
parameter space around the QCP. 

A question that comes to mind is the reason for the 
difference between our results and the ones obtained by 
MMS ||]. We claim that our approaches are valid in 
different regimes. A clear way to understand this dif- 
ference is to perform the same calculation for a Poisson 
distribution of Dirac delta potentials ||T^. In this case 
the problem is characterized by two physical parameters: 
the density of impurities p and the strength of the po- 
tential V. A universal Gaussian distribution like the one 
discussed in this paper is obtained by taking the limit of 
p ^ oo and F ^ so that u — pV^ is constant. On the 
other hand the single impurity limit studied by MMS is 
obtained by letting p — > 0. 



In summary, we have studied the problem of droplet 
formation and dynamics close to a QCP of a disordered 
itinerant Ising magnet. We find that the droplets have 
a finite Kondo temperature (associated with their dissi- 
pative quantum dynamics) that varies continuously with 
the distance from the QCP. We show that the charac- 
teristic Kondo temperature of the droplet is finite except 
in a non-universal region (dependent on microscopic de- 
tails) close to the QCP. These results are in agreement 
with our previous analysis of quantum Griffiths-McCoy 
singularities in the disordered Kondo lattice model and 
extends the validity of our results to other strongly cor- 
related systems. 
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